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ABSTRACT 

We report the detection of C0(2^1) emission from the z = 1.2 ultraluminous infrared galaxy (ULIRG) 
GOODS J123634.53+621241.3 (also known as the sub-millimeter galaxy GN26). These observations 
represent the first discovery of high-redshift CO emission using the new Combined Array for Research in 
Millimeter- Wave Astronomy (CARMA). Of aU high-redshift {z > 1) galaxies within the GOODS-North 
field, this source has the largest far-infrared (FIR) flux observed in the Spitzer 70 /im and 160 /im bands. 
The CO redshift confirms the optical identification of the source, and the bright C0(2^1) line suggests 
the presence of a large molecular gas reservoir of about 7 x 10^^ Mq. The infrared-to-CO luminosity ratio 
of L(IR)/L'(CO) = 8O±3OL0(Kkm s'^pc^)"! is slightly smaller than the average ratio found in local 
ULIRGs and high-redshift sub- millimeter galaxies. The short star-formation time scale of about 70Myr 
is consistent with a starburst associated with the merger event and is much shorter than the time scales 
for spiral galaxies and estimates made for high-redshift galaxies selected on the basis of their B ~ z and 
z — K colors. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: individual 
(GOODS J123634.53-I-621241.3) — galaxies: starburst 



1. INTRODUCTION 

Observations of molecular gas are fundamental to our 
understanding of galaxy evolution by providing measure- 
ments of the material from which stars form. The discov- 
ery of sub-millime t er galaxies (S MGs, IS mail et"al] Il997t 
iHughes et aHll998l : lBaTger et al.l ll998) enabled the ability 
to measure the molecular CO propertie s of the most lumi- 
nous i nfrared sources a t high-redshift (iFraver et al.lll998l 



1999": 'Neri et al.l l2003t iGreve et al.l l2005l : iTacconi et al.l 
200&i2008lV These observations have shown that the CO 



properties of the SMGs are similar to local ultraluminous 
infrared galaxies (ULIRGs, Ljr > 10^^ L©), and their 
number densities indicate that ULIRGs are 1000 times 
more common at z ~ 2 than locally (jChapman et al.] 
l2005f ). The importance of ULIRGs at high-redshift has 
been reinforced by dee p Spitzer mid-infra r ed (MIR) sur- 
veys at 24 ^m (e.g., ICharv et al.l |2004 iLagache et al.l 
[200I . At z ~ 1-2, 
and AGN activity 



the majority of star-formation 
occurs within dus t -rich luminous 
infrared galaxies (Le Floc'h et al.l l2005t iPapovich et al.l 
I2OO6I: ICaputi et al.ll2007fl ! 

Although the bulk of the infrared luminosity of galax- 
ies arises in the rest-frame far-infrared band (FIR, 40- 
120 /.tm), this band has been mostly unexplored at 
high-redshift. The majority of high-redshift ULIRGs 



have been identified by 24 /im and sub-mm/mm sur- 
veys. These surveys yield highly uncertain bolomet- 
ric corrections for the infrared luminosity. Deep sur- 
veys at 70 ^m and 160 /im with Spitzer enable the di- 
rect identification of luminous sources within the FIR 
band. In this Letter we report on CO observations of 
the brightest FIR source at z > 1 in GOODS-North. 
The source GOODS J123634.53-h621241.3 (also known as 
SMC GN26) was identified as b eing unusually brig ht at 
70 & 160 Aim ()Fraver et all 120061 : iHuvnh et al] 120071) and 
for having strong emission fro m polycyclic aro matic hy- 
drocarbon (PAH) molecules (Pope et al.ll2008[ ). At op- 
tical wavelengths, GN26 has a merger-like morphology 
with a disturbed core and evidence of tidal tails on 
larger scales. Although it has a large FIR flux, it is 
among the faintes t 850 jum s ource s detected to date 
(jBorvs et al.l l2003t iPope et al.l l2006t) . A cosmology of 
hm = Ho(70km s"^Mpc"^)-i = 1, 0.3, and 

r^A = 0.7 is assumed throughout this paper. 

2. OBSERVATIONS 

The C0(2^1) observations of GN26 were taken be- 
tween 2007 April 14 - June 28 with the Combined Array 
for Research in Millimeter- Wave Astronomy (CARMA) 
in the low-resolution D-configuration, resulting in a to- 
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tal of 26.4 hours of on-source data. The adopted phase 
center was the IRAC position of a(J200 0)=12^36'"34!5, 
(5(J2000)=+62°12'41" f Pope et all [2006^. The CO hne 
was observed using a digital correlator configured with two 
adjacent 13 x 31.25 MHz bands centered on 103.89274 GHz 
in the upper side-band, corresponding to CO (2^1) emis- 
sion at the repor ted optical redshift of z = 1.219 
(jCohen et al.lll996l ). The nearby quasar J1153+495 (14° 
from GN26) was observed every 15 minutes for ampli- 
tude and phase calibration, and the bright quasars 3C 273, 
3C345, and J0927-I-390 were used for passband calibra- 
tion and pointing. We estimate an uncertainty of 14% for 
the absolute calibration scale based on the observations of 
3C273 and 3C345. 

The data were reduced and calibrated using the Multi- 
channel Image Reconstruction, Image Analysis, and Dis- 
play (MIRIAD) software package (Sault et al. 1995). An 
interactive UV-plotting tool was used to visualize and flag 
spurious visibilities per baseline (105 baselines in total). 
Figure 1 shows the "dirty" (no cleaning) natural-weighted 
integrated C0(2— >1) map made by averaging over the 10 
channels showing CO emission at the phase center. No 
improvement in image quality was made using clean algo- 
rithms since the side-lobes are weaker than the noise due 
to the large number of baselines. The data are consistent 
with an unresolved point source. The weak tail extending 
eastward is not currently significant (Fig. 1). Higher res- 
olution data with higher signal-to-noise would be needed 
to measure the CO morphology of the source. 

3. RESULTS 

The total C0(2^1) line fiux of 3.45 ± 0.46 Jy km s"^ 
was derived from the integrated CO map (7.5(7, Fig. 1). 
The CO position and the limit to the CO source size of 
< 4" (33kpc) were derived from a "robust" -weighted im- 
age (ROBUST= 0). The derived position is consistent 
with previous optical and radio positions. The CO emis- 
sion peaks at the central location of the merger event. 
With the current sensitivity and resolution, there is no ev- 
idence for CO emission from the nearby comp anion galaxy 
GOO DS J123634.69-I-621243.7 at z = 1.225 (ICowie et all 
12004V which is about 3" north-east of GN26. 

Figure 2 shows the CO (2^1) spectrum at the peak of 
the CO image. The spectrum was Hanning-smoothed for 
display purposes, but a Gaussian fit of the un-smoothed 
data was used to derive the CO redshift and line width 
(Table 1). The CO redshift of z(CO) = 1.2234 ± 0.0007 is 
slightly red- ward ( 600 km s~^) of the reported optical red - 
shift oiz^ 1.219 (|Cohen et al.lll996l : ICowie et al.ll2004D . 
which is not unusual for dust-obscured ultra -l uminous sys- 
tems (iFraver et al.lll999t iGreve et alll2005f) . IWirth et all 
(|2004[ ) observed the galaxy with the Keck DEIMOS spec- 
trograph but did not report a redshift. However, visual 
inspection of the spectrum using the Team Keck Trea- 
sury Redshift Survey on-line database shows an emis- 
sion line consistent with [Oil] A3727A at a redshift of 
z — 1.224 ± 0.001, which agrees well with the measured 
CO redshift. In addition, the CO redshift is consistent 
with the redshift based on the PAH fine s (z = 1.23 ± 0.01) 
observed with Spitzer IRS (jPope et^al]|2008( ). Based on 
its sub-mm and radio flux densities (S850= 2.2 mJy; 

11 We adopt an IR definition of 8-1000 ^lra and the IRAS FIR definition 



S1.4GHz= 0.19 mJv. [Pope et all [20061 ). the continuum 
level at 3 mm is expected to be less than about 0.05 mJy. 
As expected, no continuum was detected from the co- 
addition of the lower and upper side-band line-free chan- 
nels; S(3mm)< 0.66 mJy (3ct). Since the 3 mm continuum 
level is negligible in comparison to the strong CO line, no 
continuum was subtracted from the data. 

The observed C0(2^1) line flux implies an in- 
trinsic CO fine luminosity of L'[C0(2 - 1)] = 
(6.8 ± 1.8) X 10^°Kkms~^ pc^ (see formulae 
ISolomon fc Vanden BoutI |2005[ ). which is a factor 
150-200 times larger than that for 
CO luminosity is related to the 
(including He ) by M (H^)/L^(CO) 



iTacconi et al.l 
tor a for the SMGs is 



m 
of 

our Galaxy. The 
mass of molecular gas 
a. Given that 
have found that the conversion fac- 
similar to that derived for local 
ULIRGs, we have adopted the average ULIRG value for 
GN26. Assuming a(l - 0) = 0.8 M0( Kkm s"^ pc^)"! 
and L^[CO(2-1 )1/L'[CO(1-0)] ~ 0.8 as found for 
local ULIRGs (jPownes fc SolomonI |1998[ ). we adopt 
a(2-l) = lMo(Kkm s-^pc^)-!; i.e., a(2 - 1) = 
a(l - 0)(L'[CO(2 - 1)]/L'[C0(1 - 0)])"!. For a(2 - 1) = 
1, the molecular gas mass for GN26 is M(H2) ~ 
7 X 10i°Mo. 

The baryonic gas fraction provides an estimate of the 
evolutionary state of the system h i = Mgas/(Mgas -|- 
Mstars)] fe.g.. iFraver fc Brownl 119971 ). Massive local spi- 
rals h ave gas fractions of ^ ~ 0.05-0.1 ([Young fc Scovillj 
|1991|) , while younger gas- rich systems have higher gas frac- 
tions and more evolved systems have lower gas fractions. 
To derive t he stellar mass and the ga s fraction for GN26, 
we used the iBruzual fc Charlotl (|2003 f) population synthe- 
sis models to fit the mult i- wavelength (U through z-band 
and Spitzer IRAC 3-6/^m) photometry. We assumed solar 
metallicities, a Salpeter IMF, and star-formation histories 
ranging from a single instantaneous burst to various e- 
folding timescales, including constant star-formation. We 
estimate an age of the stellar population of 1.3 Gyr and 
a total stellar mass of (2.1±0.6) x 10^^ M0 for an e-folding 
timescale of 1.0 Gyr. Neglecting HI {Mgas ^ M(H2)), the 
implied gas fraction is = 0.25 ± 0.10. 

The dynamical mass is highly uncertain given the lack 
of resolution of the CO data and the uncertain kinematics 
associated with the merger event. However, the CO size of 
0(FWHM) < 4" (diameter < 33kpc) can provide a crude 
constraint on the dynamical mass. For a wide range of 
mass distributions, the dynamical mass can be approxi- 
mated by Mdy„ « R{/W/[2sm{i)]f /G, where is the 
observed FWHM line width, G is the gravitational con- 
stant, and i is the inclination. If the intrinsic axial ratio is 
about 1, the observed optical morphology implies an incli- 
nation of i ~ 40°. For a radius of i? < 16.5 kpc, a CO line 
width of 560km s~^, and i ~ 40°, the dynamical mass is 
-^'^dyn < 7 X 10"^^ M0. This result impfies a gas- fraction 
of larger than 10% which is consistent with the estimate 
based on the stellar mass. 

We derive a FIR (42. 5-122. 5/im) luminosity of 
L(FIR) = (4.0±1.2)x 10^2^0 and a total infrared luminos- 
ity (8-1000 Atmj3 of L(IR) = (5.6 ± 1.7 ) x lO^'^Lp, using 
the 70, 160, and 850 /im flux densities (iPope et al.l 120061: 
iHuvnh et al.l [2007) and the IRS spectrum ( Pope et al] 



of 42.5— 122.5 /im throughout this paper. 
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l2008t ). A single-temperature greybody with = 41 ± 
3K and a dust-emissivity index of /3 ~ 2 fit the data. 
Similar luminosities have been previously derived based 
on local spectral e nergy di stribution (SED) templates 
(|Huvnh et al.ll2007l: IPope et a l. 2008), but a simple grey- 
body provides a better fit to the FIR data for this system. 

4. DISCUSSION 

GN26 has a disturbed, merger-like morphology (Fig. 
1) similar to local ULIRGs and many high-redshift 
ULIRGs/S MGs. The Spitzer mid-infrared IRS spec- 
trum (Pope et al. 20081 ). the infrared-to-radio flux den- 



sity ratio (iHuvnh et aLl |2007[ ). and the X-ray emission 
(| Alexander et al.ll2005f ) are all consistent with GN26 be- 
ing predominantly powered by star formation. The large 
C0(2^1) flux is consistent with this view. 

We derive a total star-formation rate of SFR~ 
95OM0yr~\ using SFR(M0yr-i) = 1.7 x 
10~^° L(IR) / Lq which assumes a Salpeter IMF (|Kennicuttl 
119981 ). The L(IR)/L'(CO) ratio provides an indication 
of the intensity of star-formation and star-formation ef- 
ficiency. Strong starbursts and ULIRGs tend to show 
high IR-to-CO luminosity ratios of > 100, while lo- 
cal spiral gal axies have lower v alues of about 10- 
50 (e.g., Sol omon fc Vanden BoutI [2005). For GN26, 
L(IR)/L'(CO) = 80 ± 30 L0(Kkm s-^pc2)-i. An av- 
erage ra tio of 350 has been re ported for high-redshift CO 
sources ([Solomon fc Vanden B out 2005). However, previ- 
ous L(I R) estimates for th e SMGs assume a temperature 
of 40 K ([Greve et al.ll2005[ ). while multi- wavelength data 
suggest a lower averag e dust temperature of about 35 K for 
the SMG popula tion (Chaoman et al."2005'; K ovacs et al.l 
l2006tlPope etai. 2006, 2008; Huynh et al. 2007), implying 
a factor of two decrease in their luminosities. If adopting 
an average temperature of 35 K for the SMGs and placing 
the data on the same infrared lui ninosity scale, we find 
that both the high-r edshift SMGs jCreve et al.ll2005^ and 
local ULIRGs (Sol omon et al.l[l997l ) have a similar average 
value of L(IR)/L'(CO) ~ 200 ± 100 Lq( K km s'^pc^)"!. 
The value for GN26 is on the low end of the range of values 
found for the more distant SMGs and local ULIRGs. 

In comparison, the recent CO observations of the 
massive ^ ~ 2 BzK galaxies (galaxies sele cted on the 
basis of their B — z and z — K colors, iDaddi et al.l 
l2007t ) suggest a similar ratio of L(IR)/L '(CO) = 60 ± 
30LQ(Kkm s"^pc2)-i ([Daddi et al.ll2008l ). Although the 
IR-to-CO luminosity ratios for GN26 and the BzKs are 
similar, the BzKs are suspected to have much longer star- 
formation time scales. For the adopted CO-to-Il2 conver- 
sion factor (a [2 — 1] ~ 1), the estimated star- formation 
(gas consumption) time scale of M(H2)/SFR = 70Myr for 
GN26 is consistent with a merger-driven starburst and the 
short time scales estimated prev iously for starburst s in lo- 
cal ULIRGs and distant SM Gs ^Solomon fc Vand en BoutI 
120051 : iTacconi et aD l2008l ). In contrast, the BzK are 
thought to have much longer star-formation time scales 
of order > 400 Myr b ased on their estimated duty cy- 
cle ([Daddi et al .^'2007^. The gas consumption time scales 
for the BzKs based on the CO data would be consis- 
tent with these results if the Galactic CO-to-Il2 conver- 
sion factor is adopted for the BzKs (e.g., a(l — 0) = 
4.8M^('Kkm s'Sc^l-^ ISanders et all ll99lD. Within 



this paradigm, BzK galaxies would be analogous to lo- 
cal spirals with enhanced SFRs and long star-formation 
time scales, while the SMGs would be analogous to local 
ULIRGs with short time scales. However, such generaliza- 
tions may be an over simplification. High-redshift galaxies 
such as GN26 and the BzKs may have a wide range of star- 
formation time scales between the intense short-lived star- 
bursts found in local ULIRGs of order 10 Myr and the long 
time scales of order 1 Gyr found for local spirals. Future 
high-resolution CO observations are needed to constrain 
the C0-to-H2 conversion factors and the associated gas 
consumption time scales before definitive conclusions may 
be drawn. 

GN26 is not representative of most SMGs studied to 
date. It has a relatively low redshift, is faint at 850 /im, 
and is warme r than other SMGs detected at z ~ 1 
([Chapman et a l. 2005). It also has the largest observed 
S(CO)/S850 ratio of any SMG to date by a factor of 5 
(in part due to the K-correction of this ratio). Previ- 
ous high-redshift CO surveys have concentrated on the 
brightest 850 /xm sources which are typically at z > 2. For 
sources at z ~ 1, the Spitzer 70 and 160 /im surveys and 
future Herschel surveys which probe the peak of the FIR 
SED may be more effective at identifying the brightest 
CO sources than the current 850 /zm survey s. The rela- 
tive i mportance of low-redshift SMGs { z 1, Wang et al] 



ig et ai.l 
il.l l2002l: 



20061 12008). "typical" SMGs {2<z< 3. lBlain et all 

Chapman et al. 20051) . and the most distant SMGs {z > 3 , 
Dunlop et al."2004': 'Wang et al.'l2007t lYounger et al.l l2007l 
2008; Dannerbauer et al. 2008) is an active topic of dis- 
cussion. I Wall et al.l ([2008) have recently argued for two 
distinct populations of SMGs: ULIRGs at z ~ 1 and more 
luminous ULIRGs at z ^--^2-3. Comparing the CO proper- 
ties of SMGs as a function of redshift and luminosity may 
help to test the hypothesis of multiple SMG populations. 

5. CONCLUDING REMARKS 

GN26 is a merger-driven starburst with a molecular gas 
mass of about 7 x IO^^A/q. Although it is the brightest 
FIR (70 /im and 160 /im) high-redshift (z > 1) source in 
GOODS-North, it has not been previously observed in CO 
given that it is among the faintest 850 /im sources. The 
CO properties of GN26 are consistent with local ULIRGs 
and the high-redshift SMGs. The derived gas fraction of 
0.25 ±0.10 based on the stellar and molecular gas mass es- 
timates suggests that about 75% of the available baryonic 
mass has already been converted into stars. Furthermore, 
at the current star-formation rate, the molecular gas will 
be depleted within < 100 Myr. The short gas consump- 
tion time scale for GN26 is in contrast to long time scales 
found for local spiral galaxies and estimates made for the 
high-redshift BzKs. 

These observations represent the first discovery of high- 
redshift CO using CARMA. Planned improvements of the 
CARMA receivers will allow for CO observations of signifi- 
cant samples of distant ULIRGs discovered by the ongoing 
Spitzer programs and future Herschel surveys. 
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Table 1 

CO(2-+l) Observational Results 



Parameter Value 

a(J2000) 12'>36"^34!51 ± 0!08 

(5(J2000) -|-62°12'41'.'2±0'.'5 

z(CO) 1.2234 ±0.0007 

Lincwidth (FWHM) 560 ± 90 km s"^ 

S'(CO)^ 3.45 ± 0.93 Jy km s'" 



L(CO) (2.7 ± 0.7) X lO^/ifcfLg 



2 
. J 

L'(CO) (6.8 ± 1.8) X lO^^hj^ K km s"^ pe 



M(H2)^ 7 X lO^QfeftfMe 



"Observed total CO(2— ►!) line flux. Uncertainty includes the 14% systematic calibration uncertainty and the 13% (7.5cr detection) random 
noise. 

''Estimated assuming a{2 — 1) = lMQ(Kkm s"-*^ pc^)"-"^. 
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Fig. 1. — The integrated CO(2^1) map of GN26 averaged over 312.5MHz (904km s^^) overlaid on an optical HST-ACS z-band image 
(grey-scale). The la error is 0.46 Jy km s~^/beam, and the contour levels are l(7x(— 2,2,3,4,5,6,7); positive contours are solid lines and negative 
contours axe dotted lines. The synthesized beam size (white ellipse) for the observations is shown in the lower left (5'.'4 x 4'.'6, PA = 80°). 
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Fig. 2. — The CO(2—>l) spectrum for GN26. The spectrum has been Hanning-smoothed. The channel increment is 31.25 MHz (90.4km s~^), 
and the 1<t error bar per smoothed channel is shown to the right. The dotted-line shows a Gaussian fit of the line profile for the unsmoothed 
data. 



